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The study of dynamics in general relativity has been hampered by a lack of coordinate independent measures of chaos. Here we 
present a variety of invariant measures for quantifying chaotic dynamics in relativity by exploiting the coordinate independence 
of fractal dimensions. We discuss how preferred choices of time naturally arise in chaotic systems and how the existence of 
invariant signals of chaos allow us to reinstate standard coordinate dependent measures. As an application, we study the 
Mixmaster universes and find it to exhibit transient soft chaos. 
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The last century has seen the emergence of three 
branches of physics which have combined to overthrow 
the Newtonian ideal of a clockwork universe. Our new 
world view encompasses an uneasy alliance of quantum 
mechanics, relativity and, to a lesser extent, chaos the- 
ory. One of the most remarkable aspects of these theories, 
in their current form, is their marked mutual incompat- 
ibility. The existence of a minimum area in the phase 
space of quantum mechanics, h, suppresses chaos, while 
the ubiquity of chaos in classical systems demands a new 
description of the semi-classical regime. The privileged 
role played by time in both quantum mechanics and chaos 
theory has led to a head-on collision with general relativ- 
ity where the choice of a time coordinate is arbitrary. 

In this letter we offer at least a partial solution to one 
of these problems, namely, the classification and quan- 
tification of chaos in general relativity. Extensive discus- 
sion of the problem can be found in Ref. so we shall 
limit our discussion to a brief review. In general relativ- 
ity both space and time are dynamical and intermixed. 
There is no such thing as the time direction. In contrast, 
chaos theory has been developed for Newtonian dynam- 
ics where time and space are absolute and the notion of a 
mechanical phase space is clear. It may appear that the 
standard tools of chaos theory can be directly applied 
to relativity in schemes such as the ADM formalism |^] 
where an explicit space-time split is made. This is not 
the case. The coordinate freedom remains in the choice 
of lapse and shift functions which describe the 3-1-1 
decomposition. The fundamentally different role played 
by time in relativity and Newtonian mechanics manifests 
itself in the coordinate, or gauge, dependence of chaotic 
measures such as Lyapunov exponents |3| . Lyapunov ex- 
ponents are the standard method for quantifying chaotic 
behaviour as they directly measure sensitive dependence 
on initial conditions. If two initially close trajectories 
separate along a given eigendirection as 



tunately, this nice picture breaks down when applied to 
general relativity. Consider the allowed coordinate trans- 
formation i — > Inr. In terms of this time variable we find 
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then A represents the Lyapunov exponent along that di- 
rection. If A > the system is said to exhibit sensi- 
tive dependence on initial conditions with a character- 
istic chaotic, or Lyapunov, timescale = 1/A. Unfor- 
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which describes the standard power-law divergence of 
trajectories found in integrable system. In particular, 
the Lyapunov exponents in this coordinate system would 
all be zero. It should be mentioned that the Lyapunov 
exponents also depend on the choice of distance measure 
in phase space and are therefore variant under spatial co- 
ordinate transformations also. From the above discussion 
it is clear that standard coordinate dependent measures 
of chaos have to be either modified, abandoned or aug- 
mented in general relativity [Q. 

In order to find methods suitable for classifying and 
quantifying chaos in general relativity we shall exploit the 
remarkable connection between chaos and fractal curves. 
The utility of these methods comes from the coordi- 
nate or diffeomorphism invariance of fractal dimensions. 
Fractal structures can be found in the phase space of 
all chaotic system. They may be uncovered by taking 
Poincare sections, plotting attractor basin boundaries or 
by finding the intersections of stable and unstable phase 
space manifolds. These methods reveal respectively, can- 
tori B , fractal basin boundaries and chaotic invariant 
sets . Actually, fractal basin boundaries are a particu- 
lar type of chaotic invariant set. The connection between 
chaos and fractals is deep. A non-chaotic, integrable sys- 
tem has sufficient isolating integrals (constants of the mo- 
tion) to fully determine the dynamics. The trajectories 
of an integrable system are restricted by these isolating 
integrals to lie on smooth manifolds in phase space with 
the topology of n-dimensional tori. In chaotic systems 
there are insufficient isolating integrals and the smooth 
tori are replaced by fractal cantori - locally the product of 
a torus and a Cantor set. The fractal dimension of a can- 
torus captures topological information j|] about a trajec- 
tory while Lyapunov exponents measure metrical proper- 
ties. The importance of such topological information as 
a qualitative sign of chaos in general relativity has been 
emphasised by Calzetta and El Hasi |^ and Dettmann 
et. al. [no|. In what follows, we discuss quantitative re- 
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lationships that relate fractal dimensions to important 
quantities such as final state sensitivity, Lyapunov expo- 
nents and chaotic entropy. 

A unified description of chaotic dynamics is possible 
in terms of chaotic invariant sets | [Tl| ] . Chaotic invariant 
sets are formed by the intersection of stable and unstable 
manifolds in phase space. The complex, and often frac- 
tal, nature of these sets follows from the result that if a 
stable and unstable manifold intersect once, they inter- 
sect and infinite number of times. Familiar examples of 
chaotic invariant sets are strange attractors and fractal 
basin boundaries. A lesser known example is the strange 
repeller responsible for chaotic scattering p^ , p^ . 

For low-dimensional systems an intriguing result has 
been found relating the fractal dimension of chaotic in- 
variant sets to the their Kolmogorov-Sinai entropy jl^ ] 
and Lyapunov exponents. The importance of these re- 
lations in general relativity follows from the diffeomor- 
phism invariance of the multifractal dimensions Dq. The 
fractal dimensions are defined by 
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where N{e) are the number of hypercubes of side length 
e needed to cover the fractal and pi is the weight assigned 
to the i^^ hypercube. The piS satisfy Y^'^J^i Pi = 1- The 
standard capacity dimension is recovered when q — Q, 
the information dimension when q — 1, the correlation 
dimension when q = 2 etc. For homogeneous fractals all 
the various dimension yield the same result. The mul- 
tifractal dimensions Dq are invariant under diffeomor- 
phisms for all q, and Z^i is additionally invariant under 
coordinate transformations that are non-invertible at a 
finite number of points |l5|] . If a strange attractor or 
repeller has a multifractal spectrum Dq with continuous 
first derivative dqDq it is said to be uniformly hyperbolic. 
Physically this means that all trajectories on the attrac- 
tor or repeller are unstable. In practice, most attractors 
and repellers are not uniformly hyperbolic and exhibit a 
phase transition at some q = qr where Dq has a discon- 
tinuous derivative. The breakdown of hyperbolicity for 
q > qT is due to the presence of stable orbits in an oth- 
erwise chaotic system. Such systems are said to exhibit 
soft chaos. 

Remarkably, it has been shown that Di equals the 
Lyapunov dimension D^ for 2-dimensional systems. The 
Lyapunov dimension is defined by 
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where Ai > > A2 are the Lyapunov exponents in each 
eigendirection and h{/j) is the metric or K-S entropy of 
the chaotic invariant set with respect to the set's natu- 
ral measure /x (see Ref. or Ref. for definitions 
of h{fi), A and /i). The metric entropy is bounded from 
above by the topological entropy H [|l6|. As its name 



suggests, the topological entropy is diffeomorphism in- 
variant while the metric entropy can be transformed to 
zero by a coordinate change. In practice, the topological 
entropy is much harder to calculate than the metric en- 
tropy, but when it can be calculated, it provides a gauge 
invariant signal of chaos in addition to the Dq^s. 

The relation Di — Dl has been rigorously established 
for certain dynamical systems and has been numeri- 
cally confirmed for many typical systems |^ . Similar re- 
lations have been conjectured to hold for n-dimensional 
systems p8| . In Hamiltonian systems conservation of 
phase space volume implies Ai — — A2. For typical 
strange attractors and repellers the metric entropy is 
given by 



hifi) = Ai - — , 
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where Td is the decay time for trajectories leaving the re- 
peller. Strange repellers are found in systems displaying 
chaotic scattering, examples of which have been explored 
in general relativity 19 2^. For dissipative systems 
with strange attractors Td — 00 and (^) reduces to the 
Kaplan- Yorke relation [|l8|. The Lyapunov dimension of 
a strange repeller is a particularly useful measure of chaos 
as it does not require a compact phase space. Most dy- 
namical systems in general relativity do not have com- 
pact phase spaces. 

With these preliminaries out of the way, we can now 
focus on the important implications eqn.(y) has for chaos 
in general relativity. The result is direct: although nei- 
ther Lyapunov exponents nor K-S entropies are coordi- 
nate invariant, their ratio is. This result allows us to 
reinstate both Lyapunov exponents and K-S entropies as 
useful chaotic measures in a fixed gauge, so long as we 
first verify that Di has a fractional value. In practice, 
it is easy to reconstruct chaotic invariant sets numeri- 
cally 1^ and measure their fractal dimension. We use 
the fractal value of Di to prove the system is chaotic. 
Of course, even for chaotic systems there are an infinite 
number of gauge choices in which the Lyapunov expo- 
nents and K-S entropy vanish. Therefore, as a physical 
choice, we exclude all coordinate choices where chaotic 
systems have Di ~ Dl = 0/0. On reflection, it is clear 
that a preferred time choice is a natural consequence of 
chaos as the K-S entropy selects a chaotic arrow of time. 
In the confines of a good gauge choice, we can compare 
quantities such as Lyapunov times to the characteristic 
timescales of other physical processes. In practice these 
are the physically important, albeit gauge variant, ques- 
tions we need to ask. 

A related method of choosing non-pathological gauges 
employs Poincare sections. Regardless of the choice of 
gauge, chaotic trajectories will appear as fractal curves 
(cantori and stochastic layers) while stable trajectories 
will appear as smooth curves (KAM tori). Pathological 
gauge choices can be defined as those which have van- 
ishing Lyapunov exponents for cantori or non-vanishing 



2 



exponents for KAM tori. At present there appears to 
be no way to quantitatively relate the fractal dimension 
of cantori to their Lyapunov exponents. Some prelim- 
inary results have recently been found relating the gap 
structure of cantori to their Lyapunov exponents p2[ , 
and hopefully new results will soon emerge to bolster 
our qualitative picture. 

For those uncomfortable with preferred time choices, 
there is another method for quantifying chaos that is en- 
tirely gauge independent. A defining feature of chaotic 
systems is their sensitive dependence on initial condi- 
tions. Generally, a range of possible outcomes can be 
assigned to a dynamical system, and each outcome has a 
basin of attraction in the space of initial conditions. For 
chaotic systems the basin boundaries are fractal, and the 
fractal dimension of the boundary provides a coordinate 
independent measure of the chaotic dynamics |]l0| , ^ . 
The quantitative importance of the fractal dimension is 
expressed in terms of the final state sensitivity f{6) [p3| . 
This quantity describes how the unavoidable uncertainty 
in specifying initial conditions gets amplified in chaotic 
systems, leading to a large final state uncertainty. The 
function f{S) is the fraction of phase space volume which 
has an uncertain outcome due to the initial conditions 
being uncertain within a hypersphere of radius 6. It can 
be shown [ESl that 



a = N ~ Df) , 
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where N is the phase space dimension and Dq is the ca- 
pacity dimension of the basin boundary. For non-chaotic 
systems a = 1 and there is no amplification of initial 
uncertainties, while for chaotic systems < a < 1 and 
marked final state sensitivities can occur. For example, 
if a = 0.1 a 50% reduction in the initial uncertainty only 
results in a 7% reduction in the final state uncertainty. 
By reversing the argument, eqn.(|^) can be used to quickly 
determine the dimension of the basin boundary . 

One application for the gauge invariant measures Dl, 
H, and f{6) advanced in this letter might be to settle 
the long running debate over the existence of chaos in 
the Mixmaster universe ||l|,|6|. At most, the Mixmaster 
will only exhibit transient chaos as its trajectories are 
asymptoically regular |24j. To show the Mixmaster is 
chaotic, we need only consider the reduction of the full 
dynamics to the discrete {u,v) map, as this map accu- 
rately describes the majority of Mixmaster trajectories 
p5[. The map is defined by 



(u„+i, w 
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1 < M„ < 2 . 



Using the standard metrical measures of chaos it is un- 
clear whether or not the (w, v) map is chaotic. For a finite 
number of iterations, n, the map has a short-time Lya- 
panov exponent given by A„ = tt^/ (6 In 2 In Umax), where 
Umax is the largest value of u visited during the n itera- 
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FIG. 1. The multifractal dimension Dq as a function of q. 



tions p5| . However, as we take the limit n — > oo to 
recover the true Lyapanov exponent we find for typical 
trajectories that Wmax ^ oo and A = 0. This regular 
asymptotic behaviour, combined with positive short-time 
Lyapanov exponents is the hallmark of chaotic scattering. 
For the (m, v) map the chaotic scattering occurs when 
1 < u < 2. In order to prove that the (m, v) map har- 
bours a strange repeller we must reconstruct its chaotic 
invariant set. For a map this corresponds to the set of 
fixed points (u„+fe,w„+fe) = for all integers k. 

By numerically generating this set we are able to find its 
spectrum of multifractal dimensions Dq. Since the (m, v) 
map is invertible, we need only consider the chaotic fu- 
ture invariant set generated by the u map. The topolog- 
ical entropy of the set is given by 

H = lim — In [No. of fixed points at order k] 

k — * OO 

= ln2. (8) 

The result H = \n2 follows from a simple counting ar- 
gument. Since H > we are assured that the (w, v) 
map is chaotic and harbours a strange repeller. It is in- 
teresting to note that the Mixmaster {u, v) map has the 
same topological entropy as the Smale Horseshoe. As 
we would expect, the compactified version of the (u, v) 
map, the Gauss map, has a considerably larger topolog- 
ical entropy oi Hq = 7r^/6/(ln2)^ jl6|. Numerically the 
topological entropy was found to converge to In 2 very 
quickly with k, differing by less than 1 part in 1000 for 
k > 10. To be on the safe side, we chose the finite ap- 
proximation k = {1..15} when calculating the multifrac- 
tal dimensions. The task of calculating Dq is made easier 
by the dense nature of the repeller for small u. We find 
the fraction of all fixed points in a given integer interval 
[1, u] to be 



F{u) = 1-2 



(9) 



so the core of the strange repeller is strongly localised 
around u = 1. For this reason, we chose to measure the 
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fractal dimension in the interval u = [1..10] as it con- 
tains 99.8% of the repeller. The above approximations 
are particularly good for evaluating multifractal dimen- 
sions with positive q as these focus on the dense regions of 
a fractal. Conversely, our truncation to finite k and u will 
lead to large errors as q —oo. The spectrum of dimen- 
sions Dq for the u map is displayed graphically in Fig. 1 . 
For reference, Dq = 0.91 ± 0.01 and Di = 0.872 ± 0.005. 
If the map is hyperbolic at q = 1 we would expect to find 
Di = 1 — l/(ATd), in accordance with eqn.(|4|). However, 
the decay of orbits from the repeller proceeds as 

N{n) = A^i exp(-n/Td) -t- N2{n + . (10) 

The power-law tail is typical for repellers punctured by 
stable periodic orbits, and leads to what is known as soft 
chaos. This indicates that the repeller is not hyperbolic 
at q = 1 and we cannot expect eqn. (Q) to hold. Since the 
repeller is not uniformly hyperbolic, it would be inter- 
esting to verify that the Dq undergo a phase transition. 
Unfortunately, the enhanced numerical uncertainties for 
q < I makes this very difficult to study. While it is dis- 
appointing that the Mixmaster is one of the rare systems 
where eqn.(|4|) does not hold, we have nonetheless been 
able to extract useful information about the dynamics 
by measuring the multifractal dimensions and topologi- 
cal entropy of the underlying strange repeller. To sum- 
marise, we have shown that the Mixmaster universe ex- 
hibits transient soft chaos due to a weak strange repeller. 
Our fractal approach has suceeded where metrical meth- 
ods have failed. 

On a grander scale, it may be that chaos theory has 
a role to play in reconciling quantum mechanics and 
general relativity in the context of quantum cosmology 
p6[ , where chaos related phenomena such as decoher- 
ence, Fokker-Planck diffusion and dynamical arrows of 
time are thought to be important. 

This work builds on my collaborations with Carl 
Dettmann, Sam Drake, Norm Frankel and Janna Levin. 
I would like to thank Albert Fathi, Robert MacKay and 
Edward Ott for answering several questions concerning 
cantori and dimensions. 
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